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I. INTRODUCTION 

The theoretical description of relaxation and decoher- 
ence processes in open quantum systems often leads to 
a non-Markovian dynamics which is determined by pro- 
nounced memory effects [l[. Strong system-environment 
couplings 2, 3], correlations and entanglement in the ini- 
tial state [J, [5|, interactions with environments at low 
temperatures and with spin baths (fj, finite reservoirs 
@, [1] , and transport processes in nano-structures 0] can 
lead to long memory times and to a failure of the Marko- 
vian approximation. 

A systematic approach to non-Markovian dynamics is 
provided by the projection operator techniques [l(J EH, 
1 1 21 ] which are extensively used in nonequilibrium ther- 
modynamics and statistical mechanics [131 ] . These tech- 
niques are based on the introduction of a certain pro- 
jection superoperator V which acts on the states of the 
total system. The superoperator V is the mathematical 
expression for the idea of the elimination of degrees of 
freedom from the complete description of the states of 
the total system: If p is the full density matrix of the 
composite system, the projection Vp represents a certain 
approximation of p which leads to a simplified effective 
description of the dynamics through a reduced set of vari- 
ables. The projection Vp is therefore referred to as the 
relevant part of the density matrix. 

With the help of the projection operator techniques 
one derives closed dynamic equations for the relevant 
part of the density matrix. The equation of motion for 
Vp can either be the Nakajima-Zwanzig equation [Tol. [TTI] . 
an integrodiffcrcntial equation with a retarded memory 
kernel, or else the time-convolutionless master equation, 
which is a time-local differential equation of first order 
involving a time-dependent generator [lj| ■ In most cases 
these equations are used as starting point for the deriva- 
tion of effective master equations through a systematic 
perturbation expansion with respect to the strength of 
the system-environment coupling. 
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In the standard approach to the dynamics of open sys- 
tems one chooses a projection superoperator which is de- 
fined by the expression Vp = ps <& po, where ps = tr eP 
represents the reduced density matrix of the open system, 
tr# denoting the trace over the environmental Hilbcrt 
space, and po is some fixed environmental state. A su- 
peroperator of this form projects the total state p onto a 
tensor product state, i. e., onto a state without any statis- 
tical correlations between system and environment. We 
emphasize that this ansatz does not imply (as is some- 
times claimed) that one completely ignores all system- 
environment correlations. It only presupposes that all 
correlations which are present in the initial state or are 
generated during the time-evolution can be treated as 
perturbations within the framework of the projection op- 
erator techniques. 

The projection onto a tensor product state is widely 
used in studies of open quantum systems. It is often 
applicable in the case of weak system-environment cou- 
plings. Usually, the perturbation expansion is restricted 
to the second order (known as Born approximation), from 
which one derives, with the help of certain further as- 
sumptions, a Markovian quantum master equations in 
Lindblad form [l5l . [lrjl [ItJ ■ In this paper the quantum 
dynamics of an open system is said to be non-Markovian 
if the time-evolution of its reduced density matrix cannot 
be described (to the desired degree of accuracy) by means 
of a closed master equation with a time-independent gen- 
erator in Lindblad form. 

A possible approach to large deviations from Marko- 
vian behavior consists in carrying out the perturbation 
expansion to higher orders in the system-environment 
coupling (several examples are discussed in Ref. [l[). 
However, this approach is often limited by the increasing 
complexity of the resulting equations of motion. More- 
over, the perturbation expansion may not converge uni- 
formly in time, such that higher orders only improve the 
quality of the approximation of the short-time behavior, 
but completely fail in the long-time limit [181 ]. 

There is however a further promising strategy: To treat 
highly non-Markovian processes in a more efficient way 
one can replace the tensor product state used in the stan- 
dard Born approximation by a certain correlated system- 
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environment state. This approach has been proposed by 
Esposito and Gaspard [lij [2(| and by Budini [2l| to de- 
rive effective master equations within second order per- 
turbation theory that describe strong non-Markovian ef- 
fects. It has been demonstrated in Ref. [13] that this idea 
can be formulated in terms of a positive projection su- 
peroperator V which projects any state onto a correlated 
system-environment state, i. e., onto a state that con- 
tains statistical correlations between certain system and 
environment states. This formulation allows an immedi- 
ate application of the projection operator techniques to 
correlated system-environment states, and to carry out 
the perturbation expansion to higher orders in a system- 
atic way. An example is discussed in Ref. [22|, where 
the master equations of second and of fourth order cor- 
responding to a correlated projection superoperator have 
been constructed. 

The application of a correlated projection superoper- 
ator implies that the relevant part V p can no longer be 
expressed in terms of the reduced density matrix alone. 
Hence, employing a correlated projection superoperator 
one enlarges the set of relevant variables to capture those 
statistical correlations that are responsible for strong 
non-Markovian effects. 

In the present paper we discuss this idea of using cor- 
related projection superoperators in the analysis of non- 
Markovian dynamics. On the basis of certain general 
physical conditions we derive in Sec. |TT] a representation 
theorem for a class of correlated projection superopera- 
tors that are appropriate for the application of the pro- 
jection operator techniques. 

A central problem of the theory of non-Markovian pro- 
cesses is the formulation of appropriate master equa- 
tions that preserve the normalization and the positiv- 
ity of the density matrix (see, e. g., the discussion in 
Refs. (HIM HE IS)- In Sec. EH we develop the general 
structure of such master equations which results form 
the application of a correlated projection superopera- 
tor. Given a superoperator that projects onto a sepa- 
rable quantum state one can construct an embedding of 
the underlying dynamics into a Lindblad dynamics on a 
suitably extended state space. Employing this embed- 
ding we derive a general class of physically acceptable 
master equations which represents a generalization of the 
Lindblad theory to the regime of highly non-Markovian 
quantum dynamics. Section IIVI contains some conclu- 
sions. 



II. CORRELATED PROJECTION 
SUPEROPERATORS 

A. General conditions 

The Hilbert spaces of the open system S and of its en- 
vironment E are denoted by TLs and TLe, respectively. 
The state space of the composite system is given by the 
tensor product TL — TLs ®TLe- States of the composite 



system are represented by density matrices p onTL satis- 
fying p > and tip = 1, where tr is the trace taken over 
the total state space. The reduced density matrix ps of 
subsystem S is given by the partial trace taken over the 
Hilbert space He, i. e. ps = trgp. Correspondingly, the 
partial trace over TLs will be denoted by trg. 

A superoperator V is a linear map O i— > VO which 
takes any operator O on the total state space TL to an 
operator VO on TL. We consider here superoperators 
with the following properties. 

1. The map V is a projection superoperator: 

V 2 =V. (1) 

It is this formal property that allows the application of 
the projection operator techniques. For an efficient per- 
formance of these techniques the projection V p should 
represent a suitable approximation of p. A natural min- 
imal requirement is therefore that for any physical state 
p the projection V p is again a physical state, i. e., a pos- 
itive operator with unit trace. This means that V is a 
positive and trace preserving map, namely p > implies 
Vp > 0, and tx{Vp) = tip. 

2. We consider projection superoperators of the fol- 
lowing general form: 

V = Is ® A, (2) 

where Is denotes the unit map acting on operators on 
Tis, and A is a linear map that takes operators on TLe 
to operators on TLe- A projection superoperator of this 
form leaves the system S unchanged and acts nontrivially 
only on the variables of the environment E. As a con- 
sequence of the positivity of V and of condition §2$ the 
map A must be As-positive, where Ns is the dimension 
of TLs (see, e. g., Ref. [27|). In the following we use the 
stronger condition that A is completely positive, because 
completely positive maps allow for a simple mathemati- 
cal characterization (see Sec. Ill B| ). 

We discuss the implications of these conditions. From 
Eqs. (JTJ) and ([2]) we get that A itself must be a projection, 
namely A 2 = A. Moreover, since V is trace-preserving, 
the map A must also be trace-preserving. Hence, we 
find that A represents a completely positive and trace- 
preserving map (CPT map, or quantum channel) which 
operates on the variables of the environment and has the 
property of a projection. A further physically reasonable 
consequence of the positivity of A and of Eq. (2]| is that V 
maps separable (classically correlated) states to separable 
states, which means that the projection does not create 
entanglement between system and environment. 

An important goal is, of course, the determination of 
the reduced density matrix ps of the open quantum sys- 
tem. Using Eq. @ and that A is trace-preserving we get 
the relation 

p s = tr B p = ti E {Vp}. (3) 

This relation connects the density matrix of the reduced 
system with the projection of a given state p of the total 
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system. It states that, in order to determine ps, we do 
not really need the full density matrix p, but only its 
projection V p. Thus, Vp contains the full information 
needed to reconstruct the reduced system's state. 

B. Representation theorem 

We derive a representation theorem for the projection 
superoperator V from the basic conditions formulated 
in Sec. Ill Al [see Eq. (fT2|) below]. Since A is supposed 
to be a CPT map one could use, of course, the Kraus- 
Stinespring representation [28|,[29| for completely positive 
maps. However, for our purposes another representation 
is much more appropriate, which will be derived now. 

We will use the following fact from linear algebra. Con- 
sider a linear operator L : V ^ V which acts on some 
Hilbert space V and has the property L 2 = L. Then 
there exist linear independent vectors \fi) and linear in- 
dependent vectors |e<) such that (fi\ej) — Sij and 

L\v)=J2\fi){ei\v) (4) 

i 

for all \v) G V. Conversely, given two linear independent 
sets and {|ej)} of vectors in V with (fi\ej) = <5y, 

then Eq. defines a linear operator with the property 
L 2 = L. Note that we neither require that the \fi) or the 
\ei) are orthogonal, nor that L is Hermitian. 

Let us apply this fact to linear maps on the Hilbert- 
Schmidt space, i. e., we take V to be the vector space of 
operators on He with the scalar product: 

{X, Y) = tr E {Xty}. 

Then we find that any linear map A can be represented 
in the form 

AX = ^2 Bi(A i7 X) =^tr B {4x}A 

i i 

with two sets {Ai} and {B{\ of linear independent oper- 
ators on He, and that the condition A 2 = A is satisfied 
if and only if (Bi, Aj) = 6^. Since A preserves the Her- 
miticity of operators, the operators A4 and Bi can be 
chosen to be Hermitian. Hence, we obtain the represen- 
tation: 

AX = ^ix E {A i X}B i , (5) 

i 

where {Ai} and {Bi} are two sets of linear independent 
Hermitian operators satisfying: 

tr E {B i A j } = 6 ij , (6) 

The condition that A is trace-preserving takes the form: 

J2{^EBi)Ai = I E . (7) 



Finally, we have to formulate the condition of the com- 
plete positivity of the map A. A given map A is com- 
pletely positive if and only if 

(I B ®A)(|V>M)>0, (8) 

where 

a 

is a maximally entangled vector in He ®He, and {|o;}} 
is an orthonormal basis for He- To evaluate condition 
([5]) we first note that 

(I E ® A)(MM) = J2 Wl ® A (i«X^D- ( 9 ) 

a/3 

On using the representation ((5]) one finds 

A(\a)(/3\) = J2W\Ai\a)Bi = £>|Af \0)B it (10) 

i i 

where Af denotes the transpose of the operator Ai with 
respect to the given basis {|a)}. Inserting Eq. (jTTJJ) into 
Eq. ((9]) we obtain 

(J B ®A)(MM) = ^^Hialifl^^l®^ 

i a (3 
i 

We conclude that a necessary and sufficient condition for 
A to be completely positive is given by the inequality 

j2 A i® B i>°- ( n ) 

i 

Employing Eqs. ([5]) and ((2|) we obtain the following 
representation for the projection superoperator V, 

Vp = Y^teE{A i p}®B i . (12) 

i 

Given observables Ai and Bi that satisfy Eqs. ©, 0, 
and (jlip . this equation defines a projection superoper- 
ator which fulfills the general conditions formulated in 
Sec. Ill Al Conversely, given a projection which fulfills the 
conditions of Sec. Ill Al there exist observables Ai and Bi 
satisfying Eqs. ©, 0, and (TTJ) such that Eq. ||T2J) holds. 
There are in general many different sets of operators Ai, 
Bi that represent a given V . If we have a particular set 
of such operators, then the operators 

j j 

represent the same projection, where u — and v = 
(vij) are real, non-singular matrices related by u T v = I. 
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C. Examples 



D. Relevant states and observables 



Within the standard approaches one considers a pro- 
jection superoperator of the form 



Vp = (tr E p) (g> p , 



(13) 



where po is any fixed environmental density matrix. Us- 
ing a projection of this form one assumes that the states 
of the total system may be approximated by certain ten- 
sor products, describing states without statistical depen- 
dencies between system and environment. The projec- 
tion (fT3|) naturally fits into the general scheme developed 
above if we take a single A = Ie and a single B = po. 
The conditions ((6]), ([7]), and (fTTj) are then satisfied and 
Eq. (TT2D obviously reduces to Eq. (TT3J) 

In the general case, a projection Vp of the form of 
Eq. (fl2f does not represent a simple product state. We 
therefore call such V correlated projection superopera- 
tors. They project onto states that contain statistical 
correlations between the system S and its environment 
E. In the following we will consider the case that one 
can find a representation of the projection with positive 
operators Ai > and Bi > 0. Equation (jlip is then triv- 
ially satisfied. Without restriction we may assume that 
the Bi are normalized to unit trace, such that condition 
([7]) reduces to the simple form J2i Ai = Ie- Under these 
conditions V projects any state p onto a state which can 
be written as a sum of tensor products of positive opera- 
tors. In the theory of entanglement (see, e. g., Ref. (30j ) 
such states are called separable or classically correlated 
[3lj . Using a projection superoperator of this form, one 
thus tries to approximate the total system's states by 
a classically correlated state. The general representation 
of Eq. (fT2|) includes the case of projection superoperators 
that project onto inseparable, entangled quantum states. 
We will not pursue here this possibility further, and re- 
strict ourselves to positive A4 and Bi in the following. 

A straightforward example for a separable projection 
superoperator is obtained through the following construc- 
tion [19|, |2l|, |22| . We take any orthogonal decomposition 
of the unit operator Ie on the state space of the environ- 
ment, i. e., a collection of ordinary projection operators 
Hi on He which satisfy 



Iljllj = SijUj, 



Then we choose: 

Ai = IL 



Bi 



tr E {UiPo}' 



It is easy to verify that with this choice the conditions 
©, ®, and pT|) are satisfied. The explicit form of the 
projection superoperator is given by 



tv E {TliPo} 



(14) 



Given a projection superoperator we define the rele- 
vant states as the states in the range of V, i. e., for which 
the relation 

Vptel = p T el 

holds. We see that these states are of the form 



Prcl 



Bi 



where the pi may be any positive matrices such that 
J2i^ T SPi = 1< Hence, the manifold of the relevant states 
is determined by the operators Bi. 

One can transfer this concept from states (Schrodinger 
picture) to observables (Heisenberg picture). A Hermi- 
tian operator O rel on the total state space is said to be a 
relevant observable if the relation 

tr{0 Icl {Vp)} = tr{0 IelP } 

holds true for all states p. This means that the expec- 
tation value of a relevant observable in any state of the 
composite system is left unchanged under the application 
of the projection superoperator V . 

Introducing the adjoint map (defined with respect 
to the Hilbert-Schmidt scalar product for operators on 
the total state space) we can reformulate this condition 



V^O xel = Orel- 

Hence, the relevant observables are those observables 
which are invariant under the application of the adjoint 
projection. From the representation (TT2]) we get: 



V^O = Y d ^E{B i O}®A i 



(15) 



from which we infer that the relevant observables must 
be of the form 



Orel = °S ® A i 



(16) 



where the % s are arbitrary observables of the subsys- 
tem S. The structure of the relevant observables is thus 
determined by the operators Ai. 



III. DYNAMICS 

A. General formulation 

The dynamics of the total system is given by 

p(t) = U(t)p(Q)U\t), (17) 
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where U(t) denotes the unitary time-evolution operator. 
Given a projection superoperator of the form of Eq. 
one introduces the dynamical variables: 

Pi {t) = tx E {A iP {t)}. (18) 

Since the Aj are positive operators, we have pi(t) > 0, 
and by use of Eq. and of the normalization tr^B.; = 1 
we find that the reduced density matrix is obtained from 

Ps(t) 

i 

The normalization condition reads: 

trspsQk) =J2 tr sPi(t) = 1. (19) 

i 

Hence, the state of the reduced system is determined by 
a certain set of unnormalized density matrices Pi(t). 
We consider initial conditions of the following form, 

p(o) = t p (q) = J2p*(v)® b *- (20) 

i 

This equation means that the initial state belongs to 
the manifold of the relevant states. As a consequence 
there is no inhomogeniety in the Nakajima-Zwanzig or 
the time-convolutionless master equation, although p(0) 
describes a correlated system-environment state. On us- 
ing Eqs. (HU), IH1), and flU we get: 

p i (t) = J2 tT E{AiU{t)p j (0)®B j tf(t)}. (21) 

3 

To be specific we assume that the index i takes on 
the values i = 1,2, ... ,n. We define a vector g whose 
components are given by the dynamical variables p^. 

Q= (pi,p2,.-.,pn)- (22) 

Equation (|21[) then defines a dynamical transformation 
of the form 

V t : g(0) ^ g(t), 

i. e., we get a one-parameter family of dynamical maps 
{V t | t > 0}, where Vq is equal to the identity map. 
For each fixed t the map Vt transforms any initial vector 
g(0) with positive components Pi(Q) > into some vec- 
tor g{t) with positive components pi(t) > 0. Of course, 
this transformation also preserves the normalization con- 
dition (|T9"j) . 

It is important to emphasize that Vt is not a quantum 
dynamical map in the conventional sense of the theory 
of open systems, simply because it is not an operation 
on the space of states of the reduced system, but rather 
a map on the space of vectors g. The transition from 
g(0) to the reduced density matrix ps(0) = J2iPi(®) * s 
connected with a loss of information on the initial cor- 
relations, such that from the mere knowledge of ps(0) 
the dynamical behavior cannot be reconstructed, in gen- 
eral. In other words, the evolution from ps(0) to ps{t) 
for t > is not a map, i. e., there is no prescription which 
assigns to each ps(0) a unique ps(t). 



B. Structure of non-Markovian master equations 

The application of the time-convolutionless projection 
operator technique leads to a closed dynamic equation 
for the relevant part of the density matrix: 

j t Vp(t)=JC t (Vp(t)), 

where K} is a linear generator which depends, in general, 
explicitly on time. If one uses a projection of the form of 
Eq. (fT2"|) this yields a system of equations of motion for 
the densities Pi{t), 

4;Pi = fci(pi,--->Pn), i = l,2,...,n, 

with an explicitly time-dependent generator K,\ for 
each i which depends linearly on the input arguments 
pi, . . . , p n . Let us suppose that these generators may be 
approximated by time-independent generators such that 
we can write 

^p i =Ki(pi,...,p n ). (23) 

The family V t then represents a semigroup of dynamical 
transformations. Of course, this semigroup assumption 
for V t does not imply that the dynamics of the reduced 
density matrix yields a semigroup. 

A typical situation in which the semigroup assump- 
tion is valid is given by a projection superoperator of the 
form of Eq. ([T3J| . where the IT project onto certain energy 
bands of the environment, describing a structured reser- 
voir. The semigroup assumption then presupposes that 
all intra- and inter-band transitions may be described by 
means of time-independent rates obtained from Fermi's 
golden rule [jj, El 111. 

Our aim is to determine the general structure of the 
generators JCi. To this end, we demand that Eq. ([2"3")) 
preserves the positivity of all components pi, i. e., given 
positive initial data Pi(0) > the corresponding solution 
satisfies piit) > for all times. A convenient way of for- 
mulating the dynamical transformation is the following. 
We introduce an auxiliary Hilbert space C™ and an or- 
thonormal basis {| i)} for this space. Then we can identify 
the vector g introduced in Eq. (|22p with a density matrix 
on the extended space Tts <8> C™ : 

= X>®KX*l- (24) 

i 

This density matrix can be regarded as a block diagonal 
matrix with blocks pi along the main diagonal. More- 
over, the reduced density matrix ps is obtained by the 
partial trace of g taken over the auxiliary space. Hence, 
the auxiliary space represents an additional degree of 
freedom which expresses the statistical correlations in- 
troduced through the projection superoperator V '. 

The dynamical transformation Vt can be viewed as 
a CPT operation on the extended space that preserves 
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the block diagonal structure. To guarantee the conserva- 
tion of positivity we therefore impose the condition that 
there exists a Lindblad generator C on the extended space 
which also preserves the block diagonal structure, i. e., 
which has the property 

£ (l>®K><*M =Y f Ki(pi,...,p n )®\i)<i\- (25) 

If such a Lindblad generator exists, the solution of 
Eq. can be written as 



'$>(0)®|i)(i| • (26) 



This implies the required conservation of positivity of the 
components Pi(t) for all times t > 0. In mathematical 
terms, Eq. (|26|) can be interpreted as an embedding of the 
non-Markovian dynamics into a Lindblad dynamics on 
the extended state space. We now demonstrate that the 
simple condition expressed by Eq. (|25|) fixes the structure 
of the generators Ki to a large extend. 

Theorem. A Lindblad generator C on the extended 
state space with the property (|25p exists if and only if 
the generators /Cj are of the special form 

£i(pi, ■ ■ ■ ,Pn) (27) 

= -i [H\ Pi ] +J2 (x$Pi*F \ {l>"x U'x-I'] 

with Hcrmitian operators H l and arbitrary system oper- 
ators R l £. 

Proof. Assume that a Lindblad generator C on the 
extended state space with the property (f2l>l) is given. As 
for any Lindblad generator we have 



C(p) = -i 



[H, g ]+J2 (RxqR\ - \ {R{Rx, gfj , (28) 



were H = Ft and the R x are operators on the extended 
state space. These operators can always be written as 
sums over tensor products: 



ij 

Rx = Y. Rl i®m\- 

ij 



Substituting these relations into Eq. and using the 
expression (|24|) for g one easily shows that 

■i k 



where 

jyik = _ i t H ik pk _ piH ik\ 



(29) 



a Pk - -piRi'Ri 



Hence, in order for condition ([25)1 to be satisfied we must 
have D lk = for all i ^ k, and, in particular, D m = ICi. 
Setting i = k in Eq. J29]) we get the form fl27]) for the 
generator ICi, where H n = H l is Hermitian. 

Suppose now that ICi is of the form of Eq. (|27|) . Then 
we define operators on the extended state space through 

h = Y.H^mi 

i 

S$ - i#®|i>01. 

where F is Hermitian. With the help of these operators 
we define a Lindblad generator by means of 

C(e) = -i [H, g] + J2 (sx'esf \ {sfs^g} 

ijX v 

It is easy to check that this Lindblad generator indeed 
has the required property <P5)) . This proves the theorem. 



C. Discussion 

Inserting Eq. (|27|) into Eq. l[23|) we obtain the master 
equation 



i [H\ P i\ 



Under the condition of the existence of an embedding 
into a Lindblad dynamics on the extended state space 
[see Eq. (f26|)], this is the general form for the equations of 
motion of the dynamical variables pi (t) , where the H 1 are 
arbitrary Hermitian operators, and the RV are arbitrary 
system operators. 



1. Examples 

Several master equations proposed recently in the lit- 
erature are of the general form of Eq. ([30]) . The simplest 
special case of this equation is obtained if we have only 
a single component p\ = p$ such that we can omit the 
indices i and j. The master equation ([50)1 then obviously 
reduces to an ordinary master equation for the reduced 
density matrix in Lindblad form. Thus, Eq. ([50)1 can 
be viewed as a generalization of the Gorini-Kossakowski- 
Sudarshan-Lindblad theorem [l5|, [l6| . 

If we choose R l x 3 = SijR\ the master equation (|30|) 
takes the form 



Jt Pi = £i(Pi). 

Hence, we get an uncoupled system of equations of mo- 
tion with a Lindblad generator Ci for each component 
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Pi. Although each component pi(t) follows its own 
Mar kovian- type dynamics, the dynamics of the reduced 
density matrix ps(t) is generally highly non-Mar kovian. 
Master equations of this uncoupled form have been de- 
rived by Budini and applied to various models for the 
dynamics of open systems in structured reservoirs [2lj |. 

In the general case an equation of the form (|30|) in- 
volves a coupling between the components pi. An ex- 
ample of such an equation has been derived in [22] from 
a specific microscopic system-environment model. The 
model describes a two-state system with ground state 
|0), excited state |1) and level distance AE, which is cou- 
pled to an environment Q. The environment consists of 
a large number of energy levels arranged in two energy 
bands of width Se (see Fig. [1}. The lower energy band 
contains N\ levels, the upper band iVjj levels. 

The total Hamiltonian of the model is taken to be 



We employ a projection superoperator of the form 



H = Hs + He + V. 



(31) 



Hs = AEo+cj- is the free system Hamiltonian, where 
a± denote the usual raising and lowering operators of 
the two-state system. The free Hamiltonian of the envi- 
ronment is given by 

Se { Se \ 

h e = Y. ^iMM + J2 ( AE + jf 712 J mm, 

ni n-2 ^ ' 

and the system-environment interaction is described by 
V = A c(ni, n 2 )(7 + \n 1 ){n 2 \ + h.c. 



The index n\ labels the levels of the lower energy band 
and n 2 the levels of the upper band. The overall strength 
of the interaction is parameterized by the constant A. 
The coupling constants c(n±,U2) are independent and 
identically distributed complex Gaussian random vari- 
ables with zero mean and unit variance. 



AE 




n, 



FIG. 1: A two-state system which is coupled to an environ- 
ment consisting of two energy bands. 



Vp = ti E {U lP } 



Pi 



in, 



1 

N~i 

f P2 



ni + {n 2 p} 
In 2 , 

N 2 



1 



-n 2 



(32) 



where Hi (II2) denotes the projection onto the lower 
(upper) environmental energy band. This is a projec- 
tion of the form given in Eq. (|14[) . It projects onto a 
state in which the environmental state Hi/Ni is corre- 
lated with the system state pi/tisPi- The second order 
of the time-convolutionless projection operator technique 
leads to the master equation (written in the interaction 
picture) 

d 72 

— pi = 7lCT + p 2 0-_ - y{cr + CT_,pi}, (33) 
^P2 = I2O-P1O+ - y{cr-0-+,p 2 }, (34) 

where 71,2 = 2tt\ 2 Ni }2 / 'Se. This is a coupled system of 
first-order differential equations for the two unnormal- 
ized density matrices pi and p2- It can be written in 
the form of Eq. ([H]) by taking H* = 0, R 11 = E 22 = 0, 
R 12 = ^/ji<j + , and R 21 = */72<7_. As has been demon- 
strated in Ref. j22l | through a comparison with numerical 
simulations of the full Schrodinger equation, this master 
equation yields an excellent approximation of the reduced 
system's dynamics. 

The second term on the right-hand side of Eq. (|33|) 
describes changes of p\ which are due to downward tran- 
sitions of the two-state system. These are only possible 
if the lower band of the environment is populated, i. e., if 
the environment is in the state H1/N1 which is correlated 
with p\ . For this reason the second term of Eq. (|33[) in- 
volves the density p\ . Correspondingly, the first term on 
the right-hand side of Eq. (|3"3"|) describes changes of p\ 
caused by excitations of the two-state system. Such exci- 
tations are only possible if the environment is in the state 
n2///V2 which is correlated with p 2 . Therefore, it is the 
density p 2 that enters the first term of Eq. (|33[) . Anal- 
ogous statements hold for Eq. (|34|) . Hence, we see that 
the transitions described by Eqs. (|3"3"|) and exactly 
conserve the total number of excitations (see Sec. MI C"2l) . 

Let us investigate the behavior of the excited state pop- 
ulation p e (t) which is defined by the expression 

Pe (t) = (l\ Pl (t)+p 2 (t)\l}. 



Equations (|35|) and 



5*W = 



lead to 

-(71 +l2)Pe(t) +7lPe(0) 

-72<l|p 2 (0)|l)+7i(0|M0)|0). (35) 



This equation exhibits a strong non-Markovian character 
because the initial data appear as inhomogeneous terms 
on its right-hand side. These terms express a pronounced 



8 



memory effect: They imply that the dynamics never for- 
gets its initial condition, i. e., that the process is governed 
by an infinite memory time. A further remarkable feature 
of Eq. ([55)) derives from the fact that the initial data on 
its right-hand side cannot be expressed solely in terms of 
the matrix elements of the reduced density matrix. This 
means that at any time t, and even in the limit t — > oo, 
the process is strongly influenced by the statistical cor- 
relations of the initial state. Hence, the influence of the 
initial correlations never dies out and is present even in 
the stationary state. 

Recently, Esposito and Gaspard [H, H3| have derived 
a master equation from a microscopic system-reservoir 
model within second order perturbation theory, which is 
also of the form of Eq. (130|) . In their derivation the index 
i plays the role of the energy e of the reservoir. The 
corresponding density matrix pi = p e describes a system 
state which is correlated with a certain reservoir state of 
energy e. If the open system represents again a two-state 
system, the master equation of Esposito and Gaspard can 
be written as 33] 

= - i[He > pe] 

+ 2^ [ 7l(£>£ )<J+Pe>(T- \(T + <J-,Pe\ 

+72(e,£ )(J-Pe'(y+ 2 l^-^+'ftl J • 

Here, 71,2 (£, s') are certain transition rates which are de- 
termined by the parameters of the microscopic model, 
and H e is the system Hamiltonian including an e- 
dependent Lamb-type energy renormalization. One eas- 
ily checks that this master equation can indeed be 
brought into the general form of Eq. ([3TJ)) by taking 

Rf' = V7i(e,e')T+ and Rf' = y/-y 2 (e, e')(T-. 

A master equation of the general form pop has been 
derived recently by Budini [321 ]. This author suggests 
introducing an extra degree of freedom U which modu- 
lates the interaction between the reduced system S and 
its environment E. Under the assumptions that E may 
be treated as a Markovian reservoir and that the dynam- 
ics of the populations decouples from the dynamics of 
the coherences of U , one arrives at the master equation 
P"TJ)) . In a certain sense the introduction of an additional 
degree of freedom U corresponds to the extended state 
space 7is €5 C n used in Sec. IIIIBI to construct the em- 
bedding into a Lindblad dynamics. An advantage of the 
present formulation is that it avoids the use of a micro- 
scopic model for the extra degree of freedom and that 
it clearly shows the basic physical condition underlying 
the master equation (j30p : This condition is the existence 
of an effective representation of the dynamics through a 
projection onto separable, classically correlated quantum 
states. 

What are the physical conditions under which the use 
of a standard projection onto a tensor product state is 
not sufficient to correctly describe the reduced system's 



dynamics for a given system-environment model? While 
this important question seems to be difficult to answer 
in full generality and certainly deserves further investi- 
gations, important hints can be obtained already by an 
investigation of the time-convolutionless perturbation ex- 
pansion P, [3] • If the two-point environmental correla- 
tion functions do not decay rapidly in time the second 
order of the expansion cannot, of course, be expected 
to give an accurate description of the dynamics. For 
instance, this situation arises for the spin star model 
discussed in Ref. [Hj], where the second-order genera- 
tor of the master equation increases linearly with time 
such that the Born-Markov approximation simply does 
not exist. The example investigated in Ref. [13] demon- 
strates that there are even cases in which the standard 
Markov condition is satisfied although the product-state 
projection completely fails if one truncates the expansion 
at any finite order. In such cases strong non-Markovian 
dynamics is induced through the behavior of higher-order 
correlation functions. Hence, one can judge the quality 
of a given projection superoperator only by an investiga- 
tion of the structure of higher orders of the expansion. 
The standard projection and the corresponding Lindblad 
equation are not reliable if higher orders lead to contri- 
butions that are not bounded in time, signifying the non- 
uniform convergence of the perturbation expansion [22| . 



2. Conservation laws and relevant observables 

For many physical models one knows certain quantities 
C which are exactly (or at least approximately) conserved 
under the time-evolution. A great advantage of the for- 
mulation by means of the master equation ([30| is that it 
allows the implementation of the corresponding conser- 
vation laws. For instance, the master equation derived 
in Ref. [l9[ reflects the conservation of the (uncoupled) 
total system energy. 

To formulate the conservation of C one chooses the 
operators A; in such a way that C is a relevant observable 
[see Sec. Ill D] . According to Eq. (jTSJ) this means that C 
can be written in the form 

i 

Then we have the exact relation 

tr{Cp(t)} = tr{CVp(t)} = J2^s{ChPi(t)}. 

i 

Hence, we can express the conservation of C on the level 
of the master equation by means of the conservation law 

j t J2 t *s{C i sPi(t)} = 0. 

i 

By use of the master equation (|50|) this yields a relation 
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between the operators H l and RV : 
l [H\& s ] (r^C^ \ {R^RtCh}) = 0. 

Thus, known conserved quantities lead to constraints on 
the choice of the operators that enter the master equa- 
tion. 

An example is given by the quantity C = <j + tj- + II2 
which counts the total number of excitations for the 
model discussed in Sec. IIII C II This quantity is exactly 
conserved under the evolution generated by the Hamil- 
tonian (f3"Tj) . Writing 

C = (7 + er_ ® III + (ff+cr- + Is) ® n 2 

we see that C is indeed a relevant observable for the pro- 
jection (f32|) . i. e., we have V^C = C. The corresponding- 
conservation law takes the form 

p e (t) + tr s p2(t) = const. 

This fact can be used to motivate the projection superop- 
erator: With the choice of Eq. ([32]) one ensures that the 
projection superoperator leaves invariant the conserved 
quantity and that the effective description respects the 
corresponding conservation law. 

IV. CONCLUSIONS 

We have discussed the theoretical description of non- 
Markovian quantum dynamics within the framework of 
the projection operator techniques. It has been shown 
that an efficient modelling of strong non-Markovian ef- 
fects is made possible through the use of correlated pro- 
jection superoperators V that take into account statisti- 
cal correlations between the open system and its environ- 
ment. We have formulated and explained some general 
physical conditions which demand, essentially, that V can 
be expressed in terms of a projective quantum channel 
that operates on the environmental variables. On the ba- 
sis of these conditions a representation theorem for cor- 
related projection operators [Eq. (fl~2j)] has been derived. 

Employing a correlated projection superoperator in- 
stead of the usual projection onto a tensor product state, 
one enlarges the set of dynamical variables from the re- 
duced density matrix ps to a collection of densities pi 
describing system states that are correlated with cer- 
tain environmental states. By means of an embedding of 
the non-Markovian dynamics into a Lindblad dynamics 
on a suitably extended state space, we have derived the 
general structure of a master equation [Eq. (|3"0|) ] which 
governs the dynamics of the pi and models strong non- 
Markovian effects, while preserving the physical condi- 
tions of normalization and positivity. A particularly im- 
portant feature of the master equation is that it is able 



to describe very long and even infinite memory times, as 
well as large correlations in the initial state. 

We have also discussed the role of known conserved 
quantities. Such quantities may be helpful to find an 
appropriate projection superoperator by demanding that 
they be relevant observables for V . Once this has been 
achieved one can express the corresponding conservation 
laws on the level of the effective description provided by 
the master equation. 

The semigroup assumption used in the derivation of 
the master equation (130|) is not really necessary. In fact, 
the present formulation can easily be extended to the case 
of an explicitly time-dependent Lindblad generator on 
the extended state space. The resulting master equation 
is then again of the form of Eq. (f3T)|) , where the operators 
H l {t) and R%(t) now depend explicitly on time. 

The result expressed by Eq. (|3TJ|) could be particu- 
larly useful for a phenomenological approach to non- 
Markovian dynamics: For arbitrary H 1 and R l £ this 
equation represents a physically acceptable master equa- 
tion because it preserves the normalization of the reduced 
density matrix and transforms positive initial compo- 
nents Pi(0) > into positive components Pi(t) > 0. We 
emphasize that the argument leading to the form (f3T)]) of 
the master equations is non-perturbative. Given a cer- 
tain projection superoperator V defining the densities pi, 
the only assumption entering the derivation is the exis- 
tence of an embedding of the dynamics into a Lindblad 
dynamics on the extended state space. 

The projection superoperators used for the description 
of non-Markovian dynamics in Sec. IIIII have a special fea- 
ture. Namely, they project any given state onto a classi- 
cally correlated state. If a non-Markovian dynamics can 
be approximated by use of such a superoperator within 
low orders of the perturbation expansion, one can con- 
clude that the true states of the total system can be rep- 
resented effectively through classically correlated states 
and that genuine quantum correlations (entanglement) 
may be treated as perturbations. On the other hand, the 
representation theorem of Sec. Ill Bl includes the case of 
projection superoperators that project onto nonsepara- 
ble (entangled) quantum states. For such superoperators 
the arguments that led to the master equation ([30)) are 
not applicable. It remains an important open problem 
to extend the formulation developed here to the case of 
nonseparable projection superoperators. Such an exten- 
sion could enable a systematic investigation of the dy- 
namical significance of genuine quantum correlations in 
non-Markovian processes. 
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